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ABSTRACT 


The VISC code is a computer automated scheme for solving the equations 
describing the viscous, radiating flow at the stagnation -point of a blunt 
body which may or may not be ablating. The code provides a basis for obtaining 
predictions of the stagnation-point heating to a body entering a planetary 
atmosphere at hyperbolic velocities. The code is written in Fortran V and is 
operational on both the UNIVAC 1108 (EXEC 8) system in use at LMSC and the 
CDC 7600 system in use at the University of California, Berkeley. This report 
gives an overview of the VISC code computational logic flow, a description of 
the input requirements and output results and comments on the practical use 
of the code. As such this report forms a "users manual" for operation of the 
VISC code. 
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1. INTRODUCTION 


The Lockheed Missile and Space Company's VISC code is a computer imple- 
mented computational scheme which provides theoretical predictions of the 
radiative and convective energy transfer to the stagnation -point of a blunt 
body during hypervelocity entry into planetary atmospheres. Although prediction 
of the convective heating rate is a direct consequence of the solution method, 
the VTSC code was developed primarily to analyze those entry conditions where 
radiative processes dominate the total energy transport. The surface heating 
rate prediction results from a solution of the coupled equations describing 
the flow of a viscous, radiating gas at the stagnation point of a blunt' body 
which may or may not be undergoing massive ablation. In developing the code, 
the objective was always to provide quantitatively meaningful heating predictions 
as opposed to identifying qualitatively significant phenomena. Accordingly, a 
serious attempt was made to provide realistic models of the dominant physical 
processes, especially as regards the spectral nature of the radiative transport 
processes . 

Historically, the VISC code is the result of an eight-year period of 

development with each subsequent effort providing an increasingly sophisticated 

accounting of the gasdynamic and transport processes. D4SC's first effort 

analyzed the relatively simple problem of an inviscid radiating flow in the 

optically thin approximation.^ This first effort used an integral method to 

solve the flow equations and the analysis applied away from the stagnat ion -point. 

The second effort continued the methodology previously developed but extended 

2 

the equations to include viscous effects. The third effort provided a major 
step forward in complexity by dropping the assumption of an optically thin 
shock layer and accounted for the spectral behavior of the continuum absorption 
and emission radiative processes.^ In addition at this point, a finite -difference 
solution to the energy equation and inclusion of surface blowing (under a simple 



chemistry model) were added. The analysis, however, was still applicable away 
from the stagnat ion -point. 

The fourth effort considered for the first time the very complex equi- 
librium chemistry which results from injection of realistic ablation products 
(e.g., carbon phenolic) into the air shock layer. ^ This realistic treatment 
of the shock layer chemistry required implementation of the FEMP program into 
the analysis. The use of FEMP, along with an even more detailed accounting 
of the spectral radiative transport, raised a computational cost barrier 
which has since prevented the application of the direct finite-difference 
analysis away from the stagnation-point. Since this fourth stage, then, the 
analysis underlying the VISC code has restricted the results to stagnation-point 
values. Contemporary with the fourth level of effort was realization- of • the 
importance of radiative transport in atomic lines and the development of a 
methodology for line transport."* 

The fifth level of effort extended the radiative transport analysis 

to include line transport. Also, with respect to the gasdynamic analysis, a 

major improvement was the replacement of the integral method of solution of, 

the momentum and diffusion equations by a finite-difference calculation. This 

placed all three conservation equations on the same footing and permitted analysis 

of "massive" blowing problems.^ Recent steps in the development of the VISC 

code continued the improvement of both the continuum and line radiative transport 

T 8 

models and its application to the Jovian entry problem. 3 Finally, the 
simplifying but nonessential gasdynamic approximations that pp, = constant and the 
construction of an inner inviscid layer near the wall were removed. Hence, 
the latest results obtained* removed all gasdynamic approximations. These 
results are from a solution which used a completely variable pp as well as a 
technique for integrating the conservation equations toward both the shock and 
_ 

We refer specifically to the results given in an unpublished paper: K. Wilson, 

H. Woodward, M. Tauber and W. Page, "Jupiter Probe Heating Rates," presented 
at the Symposium on Hypervelocity Radiating Flow Fields for Planetary Entry, 

Jet Propulsion Laboratory, Pasadena, California, January 14-15* 1972 
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wall from the interior shock layer point where the stream function passes 

9 

through zero. 

Under the current effort, the radiative transport subroutine was further 

extended to include wall emission/reflection terms. At this time we dropped 

7 

the previously employed analytic method for obtaining the flux divergence 
due to lines and replaced that method with direct numerical differencing of 
the line flux. 

The code has been used to obtain heating predictions for reentry into 
the earth's atmosphere and for probe entry into the atmosphere of Jupiter and 
Venus. Within the context of its application by a knowledgeable user, the code 
is operational on the UNIVAC 1108 (EXEC 8) system in use at LMSC and the 
CDC 7600 system in use at the University of California, Berkeley. 
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2 . MATHEMATICAL MODEL 


The analysis starts with the equations describing a viscous, radiating 
flow and effects the simplifying assumptions: 

• thin layer approximation of the y-momentum equations 

• shockwave concentric to body 

• collision transport properties appropriate to a binary mixture 

• one-dimensional radiative transport 

• restriction to an axisymmetric stagnation -point 

The resulting differential equations are available in numerous publications. 

The particular form that has been implemented in our VISC code are derived 
in detail in the Appendix to Reference 8. The three conservation equations 
all have the general form 

where T) is a density stretched normal coordinate with the values T] = 0 
at wall and 7] = 1.0 at the shockwave. 

for momentum conservation 
for species conservation 

for energy conservation 

The quantity f is the tangential velocity normalized to the value immediately 
behind the shock. The quantity c is the blown gas mass fraction. The quantity 
g is the total enthalpy* normalized to free stream total enthalpy. 

* 

For the stagnation-point problem of interest the static enthalpy equals the 
total enthalpy (to order density ratio squared) and we will simply call g the 
enthalpy. 


* Af 

z = 


z = p 2 D ^ 
^ dl] 


dg 

Z Pr dT| 
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In general the coefficients P and Q are functions of both the 
independent variable T| as well af f, c or g . Our solution method uses a 
current guess (f 1 , c 1 or g 1 ) on the dependent variable to linearize the basic 
differential equation. A new value (f*'* 1 , c i+ ‘ L or g* + ^) then is obtained by 
direct quadrature. The boundary conditions are built-in into the quadrature 
solutions. Thus this method avoids the computational difficulties associated 
with the often used technique of solving split boundary condition problems by 
guessing the function and its derivative at one boundary and forward integrating 
toward the other boundary. 

For a nonblowing problem, the boundary conditions are the known values 

of f, c and g at the wall and shock. For a blowing problem, the value of the 

T1 coordinate at which the stream function passes through zero (T] = T|*) is 

i 9 

determined from the current estimate on the velocity profile f . Chou has 
shown how the boundary conditions at the wall and shock can be combined through 
a condition which matches magnitude and slope of the dependent variables to ■ 
yield boundary conditions at T) = T|*. Then a blowing problem offers no more 
difficulty than a nonblowing case regardless of the level of blowing. 
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3. COMPUTATIONAL STRUCTURE 


A chart is given in Fig. 1 which shows the logic sequence followed in 
obtaining an iterative solution to the conservation equations. Details on the 
input/output quantities will be given in the next section. The first step 
in the VISC code computation is the determination of the properties behind a 
normal shockwave propagating at the vehicle velocity into the ambient atmosphere. 
The strong shock approximation is invoked in which the velocity behind the shock 
is set to zero. This fixes the post shock enthalpy and permits estimation of 
the post shock pressure to within (l-€) where € is the density ratio across 
the shock. A simple iteration using the FEMP subroutine to calculate the actual 
post shock density leads rapidly to a determination of all shock properties. 

The sequential solution of the conservation equations is now initiated. 

The input estimates of the velocity, blown gas mass fraction and enthalpy 
provide the data required to start the iteration process. The first step in 
the iterative cycle is to calculate collisional transport .properties. This 
is accomplished using the FEMP package set of subroutines. Along with the 
enthalpy and elemental mass fractions (as deduced from the local value of the 
blown gas mass fraction), the FEMP rountine requires a value for the pressure. 
This value is taken as a constant equal to that behind the normal shock as 
calculated above. From FEMP, the variation in temperature, density, viscosity 
and binary diffusion coefficient across the shock layer is obtained. These 
properties provide the data required to perform a solution to the momentum 
equation using the technique discussed in Section 2. Convergence of the velocity 
profile is obtained via a local iterative loop. For a fixed set of density 
and diffusion coefficient data together with the just-calculated velocity 
distribution, it is possible to determine numerically a solution to the diffusion 
equation without iteration. 
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PRINT 

OUTPUT 

DATE 


DETERMINE SURFACE 
RADIATIVE FLUX 
USING FINAL CALCULATED 
ENTHALPY DISTRIBUTION 
(FEMP AND TRANS ROUTINES) 


Fig. 1 - Computational Logic Flow 
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At this point our procedure elects to update the species number density 
distribution to reflect the change in the blown gas mass fraction profile. It 
is essential to recognize that in flows dominated by radiative transport effects, 
temperature is a crucial variable. Indeed, the temperature which is calculated 
in FEMP from the input enthalpy-mass fraction distribution must be a "reasonable" 
distribution. Having established a reasonable temperature distribution 
from the initial guess it is imperative ...to the stability of the 
solution to prevent the initial temperature profile from being 
altered to a large degree in a single step. Such a large temperature 
change often will occur if we hold the enthalpy at its initial 
value and calculate a new temperature, density and species number 
density profile using the updated elemental mass fraction distribution. 

The reason which underlies the temperature change problem is the following. 

For flows dominated by radiative transport and especially at high Reynolds 
numbers, the shear layer is nearly isothermal. However, a rapid change in 
elemental mass fraction occurs across the shear layer. Since the effective 
heat capacity of the blown gas may be quite different than that of the atmos- 
pheric gas, then the rapid change in the elemental mass fractions is accompanied 
by rapid change in the enthalpy despite the nearly isothermal character of the 
shear layer. Hence, if in the updating of the elemental mass fractions, the 
injected gas moves to a point where the enthalpy is significantly different, 
the resultant temperature is physically unrealistic and causes catastrophic 
failure of the energy equation. 

To eliminate this high sensitivity between the diffusion and energy 
equations, our approach fixes temperature at the current value obtained on the 
first pass through FEMP. Thus in the second pass through FEMP, the temperature 
is the independent thermodynamic variable and updated values of enthalpy, 
density and species number density are determined. Of course, for a nonblowing 
there is no requirement to solve the diffusion equation, and similarly, no 
requirement for a second pass through FEMP. Thus, for a nonblowing problem, 
these two stages are bypassed as indicated in Fig. 1 

Having updated the species number densities and using the current 
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temperature profile, the next step is to determine the radiative flux diver- 
gence across the shock layer. In deriving the equations describing the 
radiative transport across the shock layer, it was convenient to use the 
physical shock layer normal coordinate. Accordingly, the first task in the 
radiative transport subroutine is to determine from the T| coordinate values 
and the density distribution the physical thickness of the shock layer. 

The radiative flux divergence is then calculated at selected points and 
interpolated to yield values at all points across the shock layer. These 
flux divergence values are used along with the current velocity profile and 
Prandtl number variation to provide a new estimate for the enthalpy profile. 


A test for convergence of the enthalpy profile is made at this point 
in the computational cycle. Further comments on the level of convergence to 
be expected will be given in Section 5* If satisfactory convergence has not 
been achieved, the entire iterative procedure is reinitiated. A number of 
investigators besides the author have observed the instability which arises 
if the calculated enthalpy profile is used unaltered as the enthalpy profile 
for the subsequent iterative cycle. In the absense of a meaningful stability 
analysis (which is likely to be impossible to achieve in view of the strong 
nonlinearity of the radiative transport terms), the following expedient has 
been employed. Let g 1 be the current guess on the enthalpy profile and 

i + l 

g the resulting calculated value. Then the next guess on the enthalpy 
profile which will be used to reinitiate the iterative loop is 


/ \ i i+l 

g next = (l ' m)g + “S 

where m is a fraction less than 1.0. A practical upper limit on m appears 
to be 0.1 As m -* 0 stability will undoubtedly be achieved but at the cost of 
numerous iterations for significant progress toward convergence. Our experience 
has shown that m ^ 0.2 to 0.3 are reasonable values for blowing problems 
while m may approach 0.5 for nonblowing. 


If convergence is achieved or, as is more likely, a maximum number of 
iterative cycles have passed, then a final pass through the radiative transport 
subroutine is made. In tnis final pass the radiative flux at the wall is 
calculated using the unaltered enthalpy profile, g = g 1+ ^. This final wall 
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radiative flux value then can be compared against the value obtained from 
the enthalpy profile which existed upon the last entry to the energy equation. 
The comparison of these two flux values provides the best measure of convergence 
Preparation of data for output, which is discussed in Section 4, completes 
the computation. 
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4- INPUT AND OUTPUT 


4 . 1 Input Data 

Table 1 lists the input data required to initiate a VISC code run. 

Also given in Table 1 is the data format, FORTRAN symbol and a synoptic 
description of each input quantity. The data is segregated into card groups. 
Depending on the amount of data required, more than one punched card may be 
needed in a given group. Figures 2 and 3 are copies from a computer listing 
of a set of input data. In the following discussion we shall proceed sequentially 
through Table 1. The reader will find it instructive to examine the typical 
values shown in Figs. 2 and 3 as we move through Table 1. The centimeter- 
gram-second system of units, along with °K for temperature, has been used for 
those input quantities given in dimensional form. 

The data specified in the first card group is self-explanatory. The 
data specified in the second card group deals primarily with description of the 
free stream properties. Vehicle velocity and ambient density are self-explanatory. 
The quantity "vehicle size" is ambiguously described. This quantity is not 
necessarily the vehicle nose radius of curvature. This readily can be appre- 
ciated if we consider the case of a flat nose cylinder as our vehicle of interest 
for which the nose radius of curvature is infinite. The resolution of this 
ambiguity stems from the realization that in deriving the solution to the 
momentum equation the tangential velocity gradient is taken as unity at the 

g 

shockwave. Also, the Maslen approximation employed assumes a streamline 

Q 

curvature that is approximately that given by the shockwave curvature. 

The above assumptions require that the normalizing dimension be the radius of 
curvature of an equivalent spherical body which supports a shockwave of known 
curvature. For spherical bodies it is reasonable to assume a concentric 
shockwave and the vehicle size is indeed in spherical nose radius. For 
extremely bluff bodies (e.g. , an Apollo vehicle shape) some a priori information 
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Table 1 


VISC Code Input Requirement 


CARD GROUP 

FORMAT 

ITEM 

1 

6a6 

RUN3D 

2 

8E10 . 4 

UINF 

RINF 

XEF 

GTW 



RW 

CANV 

XBT 



TLH 

3 

5E10.4 

ALEH 

ALPHE 

ALPC 

ALPN 

ALPO 

b 

5E10.4 

CISH 

CISHE 

CISC 

CISN 

CISO 

5 

815 

NES 


IEZ 


ITGS 

10 

MDL 

LINES 

ICWT 

IDG 


DESCRIPTION 

Alphanumeric run identification 

Vehicle velocity (cm/sec) 

Ambient density (gm/cm3) 

Vehicle size (cm) 

Wall enthalpy or wall temperature , 
see ICWT, card 5 
Blowing rate (nondimens ional) 
Convergence on energy equation^ 
Weighting factor in energy equation 
iteration 

Low to high species temperature 
• switch (°K) 


Elemental 

mass 

fraction 

of 

H 

at 

wall 

Elemental 

mass 

fraction 

of 

He at wall 

Elemental 

mass 

fraction 

of 

C 

at 

wall 

Elemental 

mass 

fraction 

of 

N 

at 

wall 

Elemental 

mass 

fraction 

of 

0 

at 

wall 


Elemental 

mass 

fraction 

of 

H 

at 

shock 

Elemental 

mass 

fraction 

of 

He at 

. shock 

Elemental 

mass 

fraction 

of 

C 

at 

shock 

Elemental 

mass 

fraction 

of 

N 

at 

shock 

Elemental 

mass 

fraction 

of 

0 

at 

shock 


Number of shock layer points (max 35) 
Number of shock layer points at which 
V.q is calc, (max 25) 

Maximum number of iterations 
= 0 no radiation; = 1 with radiation 
= 0 without molecules; = 1 with 
molecules 

= 0 without lines; = 1 with lines 
= 0 wall enthalpy input; = 1 wall 
temperature 

= 0 no diagnostic 0/P; = 1 diagnostic 

o/p 
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Table 1 (Continued) 


CARD GROUP 

FORMAT 

ITEM 

DESCRIPTION 

6 

8eio . 4 

eta(i),i=i,nes 

T]-coordinate mesh points 

7 

8E10.4 

F(l),K=l,NES 

Velocity profile guess 

8 

8E10.4 

CV(l),I=l,NES 

Blown gas mass fraction profile 
(if RVW = 0, input CV(l) = 0) 

9 

8E10.4 

G(I),I=1,HES 

Enthalpy profile guess 

10 

1615 

IETZ(j),J=l,IEZ 

Subscript I of ETA(l) at which 
V-q is calculated 


- 13 - 



Vise - VISCOUS RADIATING STAGNATION POINT FLOW CODE CASE 1 

VEHICLE AND FLIGHT CONDITIONS 

VEHICLE VELOCITY » 4,000')+06 ClVSEC 

A^BJENT DENSITY » i,9000-oT GM/CM3 

VEHICLE SIZE s 1,5900*qT C f l 
AMBIENT ENTHALPY =• 0.0000 ERcS/GM 
WALL TEMPERATURE « 4,5000*03 DEG.K 

BLOWING RATE * 4,uOOO*oi NORMALIZED 


COMPOSITION BY ELEMENTAL mass FRACTION 
ATMOSPHERE F£aT SulELP 


hydrogen 

. 740 

;ooo 

HELIUM 

.260 

: ooo 

carbon 

.000 

1,000 

NITROGEN 

.000 

,000 

oxygen 

.000 

,000 


PROGRAM OPTIONS AND PARAMETERS 

5,0 PKECENT CONVERGENCE ON ENERGY EQUATIONS 

.^o weighting factor in enthalpy profile specification 

7000, LOW TO HIGH TEMPERATURE SPECIES SWITCH <OEG,K> 

27 SHOCK LAYER POINTS 

22 SHOCK LAYER POINTS FOR FLUX DIVERGENCE 

1 MAXIMUM ITERATIONS 

coupled radiative transport 

MOLECULAR ABSORPTION INCLUDED in radiative TRANSPORT 
ATOMIC LINES INCLUDED in raoiative transport 
WALL TEMPERATURE SPEefF IEO 
DIAGNOSTIC data REQLESTED 


Fig. 2 - Input Data: Vehicle and Atmospheric Parameters 
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INITIAL SHOCK LAYER PROFILES 


ETA 

F 

cv 

c 

u.ouooo 

0.1)0000 

1,00000+00 

-2 , 80000-03 

2 , 00000-01 

9, 10n00"03 

1.00000+00 

-t .20000-03 

4,00000-01 

1.90000-C2 

1.00000*00 

1 .20000-03 

5,OUOUO-01 

2, 40000-C2 

1. 00000 + 00' 

3.30000-03 

6 , 00000-01 

3.3OuO0-C2 

1,00000*00 

7. 10000-03 

*7,00000-01 

3.30000-02 

1,00000+00 

1.4Q000-02 

7,5yaoo-oi 

4 . 30U00-C2 

9,90000-01 

1.90000-02 

6,00000-01 

5.00000-02 

9,60100-01 

2 , 8QOCu'*02 

8,20000-01 

5,40000-02 

9.7UC00-01 

3 , 20000-02 

8 I 4P0U0-Q1 

5.90(100-02 

9,60000-01 

3,90000-02 

8 , 60000-01 

6,70000-02 

9,50000-01 

5.00000-02 

a, 7oooo-oi 

7. 30u00-02 

9,30000-01 

5.80000-02 

»,auouo-Ql 

*. 10000-C2 

9 , 00000-C1 

7. 60000-02 

a,9yooc-0l 

9.60000-02 

8,60000-01 

1.00000-01 

9 | OOCOO-Q1 

1,30000-Cl 

7,60000-01 

1 .50000-01 

9,05000-01 

1 . 70U00-C1 

7.0UOOO-01 

2. 0C00L--01 

9,10000-01 

2. 90U00-C1 

4,50000-01 

4 • 20000-01 

9,20000-01 

5. 30000-Cl 

o.oocoo 

8 .3Q000-O1 

9,30000-01 

5,90000-01 

0.00000 

8. «0000»01 

9 , 40000-01 

6. 5 0 0 0 0 - C 1 

0,00000 

9 . lOOCO-01 

9,50000-01 

7, lOuOO-Ol 

0,00000 

9.30000-01 

9,60000-01 

7. 70000-C-l 

0.00000 

9 , 50000-01 

9,70000-01 

R.30(j00-Ci 

0.00000 

9 » 6QOOO-U1 

9,60100-01 

F. 90000-Cl 

0,00000 

9 . 70000-01 

9,90000-Ql 

9,40000-01 

0.GUC00 

9.80000-01 

9,95000-01 

9, 70000-Cl 

0.00000 

9.90000-01 

1,00000+00 

l.fjOjOti + OO 

o.oocoo 

1 . ooooo+oo 


'l-POTNT 


Fig. 3 - Input Data: Initial Shock I^yer Profiles 
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on shockwave curvature must be known in order to obtain the effective spherical 
nose which will support that shockwave. Only in this manner will the correct 
shock layer thickness be obtained. 

The requirement, for specifying free stream enthalpy is a result of 

recent interest in applying the VTSC code to Venus entry problems. Normally 

2 

the vehicle velocity is so large that h J /2 « 1 and the free stream enthalpy 
is quite negligible. However, care must be exercised in those cases where the 
ambient gas is composed of significant amounts of polyatomic species such as 
CO^ or H^O. The FEMP routine uses an accounting system which sets the energy 
of formation of diatomic species (N^, C^, 0^, H^) equal to zero. The triatomic 
species H^O and CO^ have a large negative enthalpy in this system. Only by 
specifying that negative enthalpy in the input data is the calculation able to 
correctly soak up the amount of energy expended in dissociating these triatomic 
molecules. Neglecting the intrinsic enthalpy, the free stream enthalpy for a 
gas containing only diatomics and H^O and CO^ is 


a H 2 0 E H 2 0 + a C0 2 E C0 2 


where the energies of formation are E° ^ = 1.71 X 10^ erg/gm and E° Q = 1.33 X 10 11 
erg/gm* In this expression q and are the mass fractions of H^O and CO^. 


The wall is characterized by either a prescribed enthalpy or (normally) 
temperature according to an option given in card group 5 * If a prescribed wall 
enthalpy is input, that enthalpy must be normalized to the total free stream 
enthalpy h c 0 + 'i/ 2 - If a prescribed wall temperature is used, the value is 
input in dimensional (°K) form. 


The blowing rate is given in terms of a normalized value m = p^V^/ p^U^. 
It should be emphasized that the blowing rate is held constant during any 
given run. The blowing rate specification is further discussed in Section 5 . 


The convergence criterion t is applied by requiring 



< T 


i+1 i 

(g i+1 + J)h 

be satisfied at all shock layer points in order for convergence to be achieved. 

For nonbloving problems, convergence to a level of t = .01 (i.e., Vjo convergence) 
is an achievable goal. However, for radiative dominated problems, a 1$ 
convergence tolerance is probably too severe. Usually a few points in the shock 
layer hold-up convergence. Hence, when a Vjo convergence level is achieved 
the radiative flux at the wall is changing considerably less than 1 °jo between 
iterations. A precision of 1$ is inconsistent with the overall accuracy of 
the optical properties, transport properties and the assumption of chemical 
equilibrium which set the true error bound on the heating predictions. For 
blowing problems, convergence to a 5$ level is rarely achieved. For blowing 
problems, convergence is much more reasonably measured by comparing the wall 
flux based on the final guessed and calculated enthalpy profile. 

The weighting factor employed in stabilizing the successive iterations 
on the energy equation was discussed in the previous section. 

The FEMP routine considers two different temperature ranges. These 

two ranges are denoted as the low temperature regime and the high temperature 

regime. The low temperature regime does not contain ionized species whereas 

the high temperature regime does not contain polyatomic species. A tabulation 

of the species considered in each temperature regime is given in Table I of 

Ref. 7* Changes to that tabulation since it was published include: (l) addition 

of C , He, H^O and C0^ to the low temperature regime; (2) addition of C^, He and 
+ 

He to the high temperature regime. The VISC code user must specify a temperature 
at which the FEMP calculation is switched from the low to high regimes. 

Selection of this low to high temperature break-point must take into considera- 
tion that it should be sufficiently high that polyatomic species would not be 
present in significant amounts (< 1$) at or above that temperature nor ionized 
species present in significant amounts at or below that temperature. The 
temperature break-point will depend on the stagnation-point pressure and 
typically will vary from 5000°K at low pressures (~ 0.1 atm) to 8000°K at high 
pressures 100 atm). The solution is not critically dependent on the value 
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of the break-point temperature selected. With a little experience in a given 
problem, an appropriate value for the temperature switch will be readily 
apparent to the user. An additional requirement is that the low to high 
temperature switch be greater than the wall temperature but less than the shock 
temperature. 

The third and fourth card groups permit specification of an arbitrary 

ablation product (card group 3 ) and atmospheric (card group b) gas composition. 

A fundamental restriction on the VISC code is that the composition of ablator 

and atmosphere be comprised of H, He, C, N and 0 elements. The species included 

in the FEMP routine were selected on the basis of being the most significant 

b 7 

for a mixture comprised of these elements. 3 For both ablator or atmosphere 
the composition is given simply on an elemental mass fraction basis. The actual 
species composition of the gases entering the shock layer is not specified as 
this will be handled by FEMP. 

The fifth card group provides control options. The number of shock 
layer points is at the discretion of the user. Storage limitations require 
that the number of shock layer points be kept equal to or less than 35 * 

As was mentioned above, the radiative flux divergence is not necessarily 
calculated at every shock layer point. Storage limitations require that the 
number of shock layer points at which the flux divergence is calculated be 
equal to or less than 25* When the number of shock layer points exceeds the 
number of flux divergence points, an interpolation procedure is used to 
determine the flux divergence at the noncalculated points. 

The maximum number of iterations provides a termination point for the 
calculation when "convergence" of the enthalpy profile is not achieved. On 
the UNIVAC 1108 EXEC 8 system in operation at LMSC, the computation time for a 
full iterative cycle (i.e., solution of momentum, mass and energy equations 
including FEMP and TRANS subroutines) for a typical problem involving 30 shock 
layer points and 25 flux divergence points is approximately 1 minute. 

Computer time, then, is normally the overriding consideration in selection of 
the maximum number of iterations. 
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The remaining quantities in card group 5 are self-explanatory options. 
•With regard to the diagnostic output, there will not be a discussion of the 
diagnostic data in Section 4.2. Without extensive study of the code details, 
the diagnostic data would have no utility to the user. 

The last sequence of card groups specifies the computational grid 
(in the ^-coordinate) and initial estimates of the normalized tangential 
velocity F , blown gas mass fraction CV, and normalized enthalpy G profiles. 
Also specified are the 7]-points at which the flux divergence is evaluated, 
denoted by Q-POINT in Fig. 3* These estimates, of course, are crucial and a 
discussion of input profiles will be given in Section 5* 

In addition to the normal input listed in Table 1, the user may also 
specify the spectral emissivity of the wall to account for surface emission 
and reflection. It was felt that it would be an inconvenience to require this 
spectral emissivity data to be input for every run. In a typical application 
of the VISC code, we expect the surface properties to be fixed for a number of 

cases in which N the atmospheric and vehicle parameters are varied. Accordingly, 

\ • v 

the spectral emissivity data are specified in an internally programmed data 
statement at the beginning of TRAITS, the radiative transport subroutine. The 
data is specified in terms of the average emissivity within each of the 26 
spectral groups employed in the continuum transport evaluation. That is to say, 
the VISC code currently represents the spectral variation of surface emissivity 
by average values at 26 locations across the spectrum. The frequency values 
of the continuum groups are evident from the output discussed in the next 
subsection. For the line transport, the VISC code lifts the surface emissivity 
values appropriate to the frequency interval of the various line groups in 
exactly the same way it determines the appropriate underlying continuum optical 
depth for each line group. Currently, the VISC code is programmed with unit 
surface emissivity for each continuum spectral interval. 

4 . 2 Output Data 

Figures 4 through 13 display copies of the computer listing of the 
output data for the problem defined by the set of input data shown in Figs. 2 
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and 3* In Fig. 4, the thermodynamic properties behind the normal shock are 
shown. In addition, the density jump p /p and the Reynolds number 

Car g 

Re = p U R/^s are output. Also included on Fig. 4 is a status report on the 
s s ® 

last energy equation iteration in terms of a comparison between the guessed 
and calculated enthalpy profile. 

Figure 5 summarizes the radiative flux picture. Both the continuum 
2 

and line fluxes (in W/cm ) at the wall Q + (o) and shock Q - (l) are given 
as a function of the continuum and line group and, for convenience, the 
average spectral values (in W/cm eV) in each group were output. The final 
quantity in Fig. 5 Is the integral of the radiative flux divergence over the 
shock layer. This integral is to be compared with the sum of the continuum 
and line total fluxes (i.e., spectral summation) emergent from shock layer 
at the wall and shock. Agreement of the flux divergence integral with the 
total summation of all flux components demonstrates adequate spatial resolution 
of the flux divergence calculation across the shock layer. 

The radiative flux is calculated at the same time the flux divergence 

is calculated. This allows us to present the radiative flux distribution 

across the shock layer in terms of the flux values at each of the Q-POINT 

locations. Figure 6 shows a representative set of this radiative flux 

data. The 7] coordinate, and corresponding physical coordinate (normalized to 

total shock layer thickness) are indicated. This is followed by a listing of 

2 

the continuum and line fluxes, in W/cm , and the spectral summation. The 
spectral continuum optical depth evaluated at the midpoint of each continuum or 
line group is presented along with the flux values. 

At the termination of the radiative flux output at the last shock layer 
point, 7] = 1, the wall convective and radiative heating values are presented. 
This is shown in Fig. 7- Of course, the radiative flux value was already 
output in terms of the continuum and line components in Fig. 5 . These two 
radiative flux components are simply summed and repeated for convenience in 
the output indicated by Fig. 7* I n addition to the absolute heating values, 
a convective and radiative heat transfer coefficient, defined as the heating 

o 

rates normalized by l/2 p^U^, is presented as shown in Fig. 7 


- 20 - 



COND ITIONS BEHIND NORMA L SHCC K 

PRESSURE 3 9,B440*o6 ATmS SHOCK LAYER THICKNESS 3 1,1691*00 CM 

ENTHALPY =' n", 0014*1? epgs/ch 
TEMPERATURE 3 1,6554*04 DEg.K 

DENSITY 3 7,1800-06 GM/CM3 

DENSITY JUMP =" 9,6460- 0 2 

VISCOSITY = 5,9562-04 GM/Cm-SEC 

REYNOLD S NU *10 |R a_ 7 ,6656*05 

LAST ENTHALPY PROFILE ITERATION 

£ -y- ; p c V'” G ( OljESS ) G < CALCl 


OiOQOOO 

OtOOoOO 

i,oocoo*co 

-2.82754-03 

"2 

,82794*03 

2*00000*01 

B . 96514-03 

9, 99999-01 

-1.20394-03 

"1 

,35208*03 

4,00000*01 

1. 86543-02 

9,99992-01 

1.19622*03 

1 

,01914-03 

5,00000*01 

2,40215-02 

9,99978-01 

3. 2900i»03 

3 

, 16435*03 

6,00000-01 

3. 0l;ili~C2 

9,99934-01 

7.07877.03 

^6 

,60569*03' 

7 , 00000-01 

3, 7045o~O2 

9, 59771-01 

1 « 397latQ2 

1 

,27Q 9 2»Q2 

7,5ocuo*oi 

4 , 32 i 6i-02 

9,99516-01 

1. 8337n-02 

1 

,75257-02 

8,00000"Cl 

5.10510-02 

9.98739-01 

2. 39569-02 

2 

,44792’*02 

8 , 2OCOO-01 

5. 54 j37-02 

9,97978-01 

2.49242*02 

2 

, 8337 1”02 

8, 4Q0U0-C1 

6.12224-02 

9.96512-01 

2. 81765-02 

3 

,34427-02 

8,60000-01 

7 , JO9O4-02 

9,93266-01 

3,6398^.02 

4 

1 22553*02 

8,7p000-0l 

7,69429-02 

9,89926-01 

3.88179*02 

5 

•01527-02 

6,80000-01 

8.76592-02 

9,82916-01 

5.07054-02 

6 

, 2 J446"02 

8 , 9OCO0-01 

1,057B1- C1 

9,64630-01 

6.47155*02 

7 

,96669-02 

9,00000-01 

1.42i5B-0i 

9,04324-01 

1,01621-01 

1 

•12429-01 

9,05000-01 

1,84318-01. 

8,11344-01 . 

1.52821-01 

1 

,54127-01 

9,10000-Ct V 

3* 14375-01 

4 ,97856-01 

3. 94614-01 

3 

,92426-01' 

9,20000-01 

5,27444-01 

2,15824-05 

8.26410-01 

8 

| 4 ijo 72-01 

9,30000-01 

5,93574-01 

0,00000 

8 . 784 iq -01 

9 

,06128-01 

9 , 40000 -C 1 

6 . 54 : j 00 -Cl 

0,00000 

9,04B6n-01 

9 

.3 J359-01 

9,50000-01 

7 , l3il2-0i 

0,00000 

9 . 26192-U1 

9 

,46621-01 

9,60000-Ql 

7.71i92-ci 

0,00000 

9,46656-01 

9 

,59216-01 

9,70000-01 

8.28762-01 

o.oocoo 

9 . 58227-01 

9 

,69564-01 

9,80000-01 

B • 86^83-ci 

0,00000 

9,68507-01 

9 

, 7'l382"0l 

9,90000-Ot 

9,43174-01 

0.00000 

9,78609-01 

9 

,80688*01 

9,95000-Cl 

9 , 71624 - qi ” 

0,00000 

9,88680*01 

9 

,92173-01 

1 , 00000*00 

1 , ooooo+oo 

0,00000 

9.98764-01 

1 

,ooooc*oo 


Fig. 4 - Output Data: Normal Shock Conditions and Status of Enthalpy Iteration 
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Pig. 5 - Output Data: Radiative Flux at Wall and Shock 
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Fig. 6 - Output Data: Typical Radiative Flux Data 

at Selected Shock Layer Points 
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Fig. 7 - Output Data: Convective and Radiative Fluxes 



A summary of the final values of the gasdynamic quantities — stream 
function, tangential velocity, enthalpy, ablation product mass fraction, flux divergence 
(w/cm^) and a repeat of the Q- POINT location — is output next. A representative 
set is shown as Fig. 8. A quick comparison between the initial and final 
profiles is thus available by examining the data of Figs. 3 and 8 . It will 
be observed that the enthalpy tabulated in Fig, 8 is the weighted average of 
the final guessed and calculated values shown in Fig. 4 . Thus this enthalpy 
profile represents the input to the next iteration. In principal, the output 
shown in Fig. 8 can be used directly to initiate a new set of iterations. 

A summary of the thermodynamic and collisional transport properties 
consistent with the gasdynamic data given in Fig. 8 is the next data output. 

A representative set is shown as Fig. 9 - 

The next series of output presents the molecular and atomic species 
number density distributions across the shock layer. A representative set 
for the molecular species is shown as Fig. 10 and likewise for the atomic 
species as Fig. 11. 

The final computation which the VISC code performs is an evaluation of 
the radiative flux at the wall using the unaltered calculated enthalpy profile 
(as opposed to the profile weighted between guessed and calculated values) 
from the last iteration. The representative values of output which describes 
this final calculation is shown as Fig. 12. In Fig. 12, data from the FEMP 
routine is displayed with the temperature variable being of maximum interest. 

The' final output gives the continuum and line flux which may be compared with 
the data shown on Fig. 5 to assess the level of convergence of the overall 
radiative flux determination. 
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Fig. 8 - Output Data: Shock Layer Gasdynamic Variables 
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Fig. 9 - Output Data: Shock Layer Thermodynamic 

and Transport Variables 
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Fig. 10 - Output Data: Shock Layer Polyatomic and Diatomic Species 
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Fig. 11 - Output Data: Shock Layer Atomic Species 



LS l MG L*3T C * l C U _ A*5 2 ? ,▼ H i. L oY JKOPJuE 


ETA 


T 

R*TJ 

MU 

R-* 


c=>n 

5k \ 

hr 

,-IH 

0 



tCEG.K) 

(gv:?;3) 

<CH/*ri-SEC) 

(GM2/CH4- 


<F.3GS/GN- 

(EHG&'C"- 





— 

- 

— 

— 


- 

.OiC. K>.__ 

S^C-il ‘vGiKl 

- “ 


— 

1 20 

1 

8.5QCO*-03 

8 , ^44 4 *?4 

1,2358-03 

1,00^8-:* 

1 

, 0?«99+(j8 

2 1 9o j6 + t]5 

4 , 457 : 1-01 

3 , • 1 5 5 1 + 1 1 

-2.G773-.i3 

• 20 

1 

4.60*0*03 

7 , 5 7 i 6 • ^ 4' 

l.?634-o3 

7 . 59 10- j7 

1 

,24;j!'+u H 

3 , 4^;;5+.j5 

*,5 4 25- )i 

2, >427 + 1-1 

-1,3521-03 

. 4 C 

l 

4 . 7455*33 

6. 99,4-. ’4 

. 1 , 3 , 334-03 

9.u?«-:7 

i 

,5i3 1*08 

4 , 16‘1+ 

4 , 7376-01 

2. 76‘5 7 + i - 

1, 0091-03 

.Sfl 

i__ 

4,8469+33 

4 . 50^5-? i 

1.353?-'j3 

j , 6 7 : 2 - : 7 

1 

, 7 ^ , l i 7 + q «i 

4 , 6 3 j 7 + j '3 

4.V332-UI 

2.'*2<i7 + .'X 

3.1643-13 

» 6 0 

1 

4*9962+03 

5 . 

1.37dl-j3 

8,00O7= : 7 


,0'j?3*08 

'1,20.52* j5 

.3,3023-01 

2.4213 + :. L 

6 , 60!?7-.j5 


1 

5,2664*33 

4. 3370»'J4 

1 • 4^4j-o3 

7.0332--J 7 

~\ 

c. 

,2-43 + Qrt 

^,35 50*3^ 


2. 1233+'. 1 

1,2709-12 

.75 

1 

5,4354*33 

4 , 23 -3 •04 

1.5U6-03 

6,40^6*; 7 

n 

•24?°+0H 

3 # l ;L 76 + C^ 

.3,6207-05. 

1.9218*61 

1,7526-«J2 

.80 

1 

5,6972*33 

3 . 52*5-o< 

1.6ij02”J3 

5,64j0-:7 

2 

» 0364+0* 

4.30-.6-0f» 

7 , 6393-01 

;t* 6/39 + - \ 

2.4470-02 

.82 

1 

5,8432*33 

3. t982-34 

1 • 632i"U3 

‘3.27 ' 7 * : 7 

1 

» 8522+1)8 

3,72 V7*,p 

3,2184-0! 

1 . i.O 7 7 + c 1 . 

2.6337-ji| 

,84 

1 

6,0773*33 

2* v 026-0 4 

1.7?8?-.j3 

4,843J s -:7 

1 

,4946+08 

2,:'6 5«* -.5 

9,0174-01 

;l, n?7+v:i. 

3,3443-^2 

~ 786 ~ 

1 

"o i 6*57*33" 

2* 21^5» ?4 

i.9 : 328-03 

4 * 2 1 J 6 3 7 


1 672i 4 ‘07 

1,45 53 * :,5 

1 , f.l : J29+'J0 

l.?305*Ci 

4.2755-0?. 

,87 

2 

7 ,7499*33 

1.7370*04 

i 2,C53i-J3 

3.36^2=:7 

T 

,5PQO+o7 

8 , 4 1 ,4’** T4 

4,782 7 -Ql 

1 ‘. 1221 +r-i. 

5 , •'•153-02 

,88 

2 

9,6463*03 

1 . 2775«04 

1 2.0{6q*u3 

2 . 5 7 4 c ; 7 

4 

tQ"i32+07 

9,57 73 + j4 

0,4^51-01 

r 233 +r 1. 

0 . 2045-02 

’.89 

2 

i, 1153*34 

' 9. 4Ci 8 * 0 ^ 

1 I. 7068-03 

1 , 6o^^ w ;'7 

0 

1 475^+c 7 

1 « 3l4V*c5 

\ * 4154-0! 

». 7 :5«V3+ f -iJ 

/,966 7-.j>; 

• 9 0 

2 

1,1793*34 

6,27*7-0* 

! 1,4734-03 

9,251X508 

9 

,5754+07 

1,7b 16+05 

11 

6, 1.753 + OJ 

1 , 1 2 4 3 - i 1 1 

.’1 

2 

1,1924*34 

4.3445*?5 

1 1.3*520-03 

5,a7^ /*0 8 

1 

,2782+08 

2,1b )5*0- 

(i tC 34^-0! 

4 • 3 2'.) 2+ C O 

1,5413-01 

«9l 

2 

1,4222*04 


T9 .0219-04 

i. 52>5 sj6 

3 

, 99;j8 + o6 “ 

"2. 9q27*'5 


2-.)047 + i:(J 

3, 92 4 3-01 

.’2 

2 

1,5523*34 

8. t 676- '* 

i 6 . 8 75 q- j4 

3 * 6 1 - 2 = 9 

9 

1 13.1.7 + 0 6 

3. ?.o 24 + ^5 

lfV2d2+U0 

1* ’372 + . 1.1 

A , 4007-ni 

.93 

2 

1,5979*34 

7. 72?6- 16 

6,4342-04 

4 ♦ 9843« 

9 

,8657+oH 

3,22 

1 . 9 7^K + ,jO 

l.-J2*i6+nj 

V , 0*13-01 

•,94- 

2~ 

IV 6119*04 

7,9805^-^ 

0,32^3-04 

4 . 0q3 ' 9 

1 

, ?n76 + c9 

3,2i-t5* ;5 

9833+00 

i**.»i96+f. n 

V. 3n36-.»l 

.95 

2 

1,6225*04 

7 , 48 * 4 *: * 

6 . 335Q-04 

4,6674-;9 

1 

,0233+UV 

3, 2042*05 

1,99 ;H*rJ0 

1 . 0126+:; U 

9,4*62-1.11 

.94 

2 

1 » 63C 7*34 

7.4074*^6 

0.1444-.J4 

4 , 567 4* 9 


i0349+oV 

3,19',O+0 s 

x .9 l, 57+ in 

1. >072+ r H 

5W-0i 

".97- 

2 

it 6 37 2+ 3~4 7Y 3 4 A 3 • 6 ‘ 

" 6,1o85-u4 

4 • 4 $ 6. 7 - r 9 

X 

.0441+09' 

5, 1R’5?.*C5 

l,9''94+mj ’ 

1.0029 + f.u 

9 ,6^56+0 1~ 

.98 

2 

1,6427*04 

7 , ?976-n6 

6,0624 + 1)4 

4»42*i°0 9 

i 

•C515+0V 

3, 1 6 J <t * ,") 5 

2.002X+JO 

9, *> 9 2 9 — ,r > 1 

9 , 7 A 3 8 - ' 1 1 . 

.99 

2 

1,64J7*04 

7,9976-06 

6,0624+04 

8,42U-j9 

X 

,Q5i5+0V 

J, ln.i4-05 

2i()O2??+00 


9,b66 y -*U 

“ .99 2 

1,6505*04' 

7.2247-!6 

5, 9065-04 

8,332J«:9 

1 

•CS21*09 

3,1 //•.„* 35 

2,i;.(jb9 + dQ 

9 * 9406-X « 

9,9217-01 ' 

1.00 

2 

1 1 65Q6+Q4 

7 9 2247*06 

5,9965+04 

4. 3323-0’ 

i 

•0621+07 

3,1746*05 

2,0039+00 

9.9406*M 

1 , 0000+00 


CELT A* 1 , 16936 l J5'*O0CM 

QONTJNUUH FUJX •_ _2,43S7*04 

.UINE'>u)»t ~ 6~ 220 0*0 3 


Fig. 12 - Output Data: Final Temperature /Enthalpy 

Profile and Resultant Surface Flux 



5. PRACTICAL ADVICE TO THE USER 


Set forth in this section are some practical guides to using the VISC 
code as deduced by the author as a result of many attempts (successful and 
unsuccessful) at obtaining heating predictions. 

The fundamental problem is prescribing a set of initial profiles 
(velocity, blown gas mass fraction, enthalpy) in a nonphysical coordinate 
which are sufficiently close to reality to get underway with a series of 
iterations. For a nonblowing problem this is not difficult at all. The 
boundary conditions on the velocity profile are f = 0 at the wall and f = 1.0 
at the shock. A linear rise from the wall value to a level f = 0.4 at the 
edge of the boundary layer is reasonable. The edge of the boundary can be 
estimated as T| (edge) ~ 100/(Re s ) 2 . From the edge of the boundary layer to 
the shock, another linear profile is reasonable. Thus an input velocity 
profile might look like that sketched below. 



V 
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The boundary conditions on the enthalpy profile are g = S wall at the 
vail and g = 1.0 at the shock. The wall enthalpy can be prescribed directly 
or, if the usual vail temperature option is being used, the vail enthalpy 
can be estimated. The first iteration will correct any error in the vail 
enthalpy estimate. Again, a linear rise across the boundary layer region to 
the edge value is an adequate initial guess. The appropriate edge value will 
depend on the degree of radiative cooling. However, the edge value will 
typically be around S 0(i g e ~ 0.8 to 0.9* It is normally best to err on the 
low side. A linear variation from boundary layer edge value to unity (the 
value at the shock) should be adequate to complete the initial enthalpy 
profile guess. Thus an input enthalpy profile might look like that sketched 
below. 



For a blowing problem, the selection of reasonable initial profiles 

is much more difficult. The user must first arrive at a reasonable bloving 

rate keeping in mind that the bloving rate will be fixed during a given 

computer n run. " In previous years ve were content to consider the bloving 

6 T 

rate as an independent parameter. 3 In our recent studies of the Jovian 

8 

entry heating problem, a strong attempt has been made to determine the bloving 
rate based on the vail radiative flux and an effective energy of ablation. 

A conservative estimate to the wall flux will be the value deduced from a 
nonbloving calculation. A correlation of the reduction in the radiative flux 
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due to blowing (e.g. , the correlation presented by Page 10 ) can be used to 
obtain a much closer estimate of the ’’final" blowing rate. 

Having fixed upon a blowing rate, the user next must decide on an 
approximate location of the extent of the blown gas layer, i.e., the value of 
the T) coordinate where the stream function is zero, T] = T|*. The shear 
layer with its rapid change in velocity, blown gas mass fraction and enthalpy 
will be centered (approximately) at this blown gas-atmospheric gas interface. 

The author can offer no general rule for determining the value of T|*. It 
will be necessary for the user to experiment with the problem at hand to 
obtain a feeling for the interface location. A single pass through the entire 
iterative cycle or even a single iterative solution to the momentum equation 
will provide a good second estimate for the value of T]*. This is true even 
if the initial guess for the profiles is so poor as to cause the VISC to 
undergo an error termination in the radiative transfer or energy equation 
routines. The user should call for the diagnostic output and determine from 
the output the value of the T] coordinate where the stream function (identified 
as PSl) passes through zero. 

Assume, then, that the user is unable to arrive at an a priori estimate 

for the interface location 1 ] = T|*. A practical approach might be to set 

T|* = 0.5 and select initial velocity, blown gas mass fraction and enthalpy 

profiles around that value. Typically, the velocity function f increases 

linearly from the wall value f = 0 to the level f = 0.1 at the inner edge 

of the shear layer. As in the nonblowing problem, the velocity at the outer 

edge of the shear layer is f = 0.4 and the velocity increases in a nearly 

linear fashion to its shock value f = 1.0. The extent of the shear layer is 

(as for the nonblowing case) on the order of 100/(Re ) s . A simple straight 

s 

line interpolation between the value at the inner and outer edge completes 
the initial velocity profile guess. A typical guess for the velocity profile 
based on T]* = 0. 5 is shown in Fig. 13 . 

The blown gas mass fraction can be taken as a constant value of unity 
(the wall value) out to the inner edge of the shear layer. The blown gas 
mass fraction can be assumed to drop in a linear manner to a value of zero 
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NORMALIZED FUNCTION (VELOCITY, MASS FRACTION, ENTHALPY) 



DENSITY STRETCHED NORMAL COORDINATE, tj 


Fig. 13 - Representative Profiles of Velocity, 
Mass Fraction and Enthalpy for a 
Massive Blowing Problem 
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(the shock value) as the outer edge of the shear layer. Such a guess for the 
blown gas mass fraction is shown in Fig. 13 . 

No general rules can be given for selection of the initial enthalpy 
profile. A crude estimate is to assume that the enthalpy is constant at the 
wall value out to the inner edge of the shear layer and then rises linearly 
across the shear layer (as in the nonblowing case described above) to a value 
on the order of 0.8 to 0.9 at the outer edge of the shear layer. Again, a 
linear variation for outer shear layer edge to the shock value of unity com- 
pletes the enthalpy guess. The enthalpy profile estimate using this approach 
is shown in Fig. 14. Under conditions of massive blowing and strong radiative 
coupling, the blown gas will experience a substantial rise in enthalpy before 
it reaches the interface. If the user initiates the solution with a constant 
enthalpy in the blown gas layer, it will probably be necessary to perform a 
number of single iteration runs observing the inward movement of the region 
heated blown gas before a reasonably consistent enthalpy profile would be 
achieved. 

The mesh selection remains to be discussed. The VISC code does not 

provide an automated routine for selection of points in the T) coordinate 

system. The user should provide for the location of a few mesh points within the 

shear layer region. It is not at all necessary nor practical to provide a 

fine mesh inside the shear layer when the Reynolds number becomes large, say 

Re ~ 10^. The user must judiciously select the remaining points to provide 
s 

resolution of the region near the shockwave where the radiative cooling 
produces a substantial drop in enthalpy and in other regions (particularly 
in the blown gas layer) where the radiative flux divergence changes rapidly. 
Normally the spatial resolution around the interface is dominated by require- 
ments to capture the peak in absorption of the intense incoming radiative flux. 

The selection of the points at which the flux divergence is calculated 
(the Q-POINT distribution shown in Fig. 3 ) reflects the comments given above 
concerning the need to obtain an accurate accounting of the rapid change in 
flux divergence which occurs in various regions near the edge and interior to 
the blown gas layer. An T| point and its designation as a Q-POINT is required 
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very near the shock, e.g., at T] = .99 or T| = *995* This requirement stems 

from the fact that the flux divergence due to atomic lines changes by at least 

one order as the calculation moves from a point exactly on the (discontinuous) 

shockwave to a point which is a minute physical distance inside the shock 

7 

layer. To avoid the requirement of a detailed resolution of the line flux 
divergence near the shock; the code uses the flux divergence from the point 
nearest the shock to provide an extrapolated value at the shock. For this 
reason, a point reasonably close to the shock corresponding to a physical distance 
of about 1$ of the shock layer thickness) is required. 

At what point does the heating prediction obtained by the VISC code 

constitute a "solution" to a particular entry problem? This question is at 

the heart of the convergence difficulties. For nonblowing problems the author 

has normally set a convergence level of 5$ with the realization that the 

precision in the radiative and convective flux calculation will be at least 

that good. For a blowing problem the author has worked toward a goal of a 

10$ level of precision. By this is meant the following. Let the radiative 

til 

flux at some cycle be the i value. Then, if the computation were continued 
an endless number of iterations, the asymptotically approached radiative heating 
value would be within 10$ of the i value. The degree of precision is measured 
by comparing the radiative flux. from the guessed and calculated enthalpy. 

The user is warned, however, that in order for the above measured precision to 
be valid the enthalpy profile must be converged to about 10$ at most points 
across the shock layer. The exception is the few points in the shear layer 
region where large spatial variations occur: In the final analysis it is the 

responsibility of the user to interpret the results of the VISC code in an 
intelligent way. 
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